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former is being used in calculations of electroweak matrix elements, while the latter have been used 
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I. INTRODUCTION 

Improved staggered fermions are an attractive choice 
for the numerical study of QCD, and are being used for 
a variety of calculations relevant to phenomenology. For 
calculations of electroweak matrix elements, such as our 
ongoing calculation of Bk [H-Q , one needs to match con- 
tinuum operators in the effective Hamiltonian onto corre- 
sponding lattice-regularized operators. Here we calculate 
such matching factors for fermion bilinears composed of 
improved staggered fermions with various gluon actions. 
We work at one-loop level in perturbation theory. 

The motivation for this work is three- fold. First, 
the results are a step on the way to the calculation of 
matching factors for four-fermion operators, such as that 
needed for Bk, results for which will be presented in an 
upcoming work Second, our results allow us to com- 
pare the efficacy of improvements to fermion and gauge 
actions at reducing matching factors. Third, some of 
our results can be compared to ongoing calculations of 
matching factors @ using non-perturbative renormaliza- 
tion (NPR) ^7^. We can also check our result for the mass- 
renormalization for asqtad fermions with that obtained 
(as a byproduct of a two- loop calculation) in Ref. Q. 

Two major problems with unimproved staggered 
fermions are large taste-symmetry breaking and large 
perturbative corrections to matching factors. Previous 
work has shown that both problems are alleviated by 
smearing the gauge links to which the fermions couple. In 
particular, it turns out that HYP smearing[9]^ is most ef- 
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fective at reducing one-loop perturbative corrections [111 ] , 
and also in reducing the taste symmetry-breaking in the 
pion spectrum [l^, U^l • In light of this we are using such 
smearing for valence quarks in our ongoing calculations of 
matrix elements. These calculations make use, however, 
of the MILC configurations [3], which use a Symanzik- 
improved gauge action. Thus we have undertaken the 
extension of the calculation of matching factors to the 
improved gluon action. We have done so using both the 
HYP-smeared action but also using the asqtad action. 
The latter gives further information on the comparison 
between smearing methods. 

The generalization to an improved gluon action is non- 
trivial. The gluon propagator is diagonal with the Wilson 
gauge action (in the Feynman gauge) , but becomes a full 
4x4 matrix for an improved gauge action. Thus various 
simplifications that are possible with the Wilson gauge 
action do not occur with the improved gauge action. 

A calculation along similar lines has been done pre- 
viously in Ref. [l^. The authors consider the asq- 
tad fermion action and Symanzik-improved glue, but 
use staggered fermion operators containing unsmeared 
("thin") links. This is in contrast to the operators which 
we use, in which all links are smeared. 

The paper is organized as follows. In Sec.|ITl we explain 
our notation and conventions for actions and operators. 
In Sec. mil we describe the renormalization of bilinear 
operators on the lattice. In Sec. IIV[ we explain the pro- 
cedure of matching between the continuum and lattice 
operators. In Sec. |Vl we close with a discussion of our 
numerical results. We relegate technical details to three 
appendices. Appendix VK\ discusses the gluon propagator 



links introduced in Ref. [To|l . We refer to these simply as HYP 
links in the following. The two smearings differ at higher-order 
and non-perturbatively. 
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for improved actions, App. |B] gives results for the renor- 
malization of HYP fermions, and App. [C] describes the 
resuhs for asqtad fermions. 

A prehminary account of this work has appeared in 
Ref. dl]. 



II. ACTIONS, FEYNMAN RULES AND 
OPERATORS 



If we consider only on-shell improvement, then one op- 
erator is redundant [l8l |. and we adopt henceforth the 
convention of setting C3 = 0. 

The Wilson gauge action corresponds to the choices 
cq — 1, ci — C2 — 0. As shown by Liischer and Weisz, 
tree-level on-shell Symanzik-improvement of the pure- 
gluon theory is obtained if the improvement coefficients 
takes the values [13, fisl 



A g eneral form for the 0(a^ )-improved gluon action 
is jl7l | (using the labeling convention of Rcf. [l8| ) 



So 



9l 



co^iReTr(l-[/pi) 



pi 



ci ^I^eTr(l - U,t) 
C2 iReTr(l - U^^) 



C3 



Y ^ReTr(l - f/eh) 



ch 



(1) 



Here pi, rt, pg, and ch denote the shape of the Wilson 
loops — plaquette, rectangle, parallelogram and chair, re- 
spectively. The overall normalization of the coefficients 
is such that 



Co -I- 8ci + 8c2 -f I6C3 = 1 , 



(2) 



Co 



ci 
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C2 =0. 



(3) 



The MILC collaboration use a one-loop improved action 
determined in Refs. [l^ [l^l- In a perturbative calcula- 
tion, however, the one-loop corrections to the improve- 
ment coefficients enter at two-loop order in a calculation 
of bilinear matching coefficients. Thus, for our one-loop 
calculation of matching coefficients, the consistent choice 
is to use the tree-level coefficients ([3]) when determining 
the gluon propagator. 

The propagator for the improved gluon action is well 
known. We have found a relatively simple form for this 
propagator, which is given in Appendix \K\ 

The staggered fermion actions that we consider in this 
paper are the asqtad and HYP actions. The former is 



'^asqtad 



\[W^(n)x{n + A) - Wlin - ji)x{n - /t)] , 

^[?/m(")x(" + A) - Ul{n - fi)x{n - fi)] , 

-^3 [U {n, n + 3fi)x{n + 3fi) - U {n, n - 3fi)x{n - 3/i)] , 



(4) 

(5) 
(6) 

(7) 



r 



where n = (ni, ^2, ^^3, ^^4) labels lattice sites, rjf^{n) = 
(— 1)"!^ is the usual staggered phase, U^{n) is 

the original, "thin" link, and m is the quark mass in 
MILC's convention. Here and in the following we set 
the lattice spacing to unity, except where clarity dictates 
otherwise. W^{n) is a smeared link constructed using 
the Fat7 blocking transformation [2l|, [l^] combined with 
Lepage's prescription 23 1 and tadpole improvement [23 |. 
U{n, n±3/i) are products of 3 thin links in the /i direction. 



U{n,n + ifi) = C/,,(n)t/^(n + /i)[/^(n + 2A) 
C/(n,n-3A) = Ul{n-ii)Ul{n-2fi)Ul{n- 



3A) 



and appear in the Naik term. Finally, wq is the tad- 
pole improvement factor, which we determine here as the 
fourth-root of the average plaquette. This action is tree- 
level 0{a^) improved. 

The HYP action is simply the unimproved staggered 
action using HYP-smearcd links: 



5hyp = Y.^{n)[Yri,{nWl 



X{n) , (8) 



V"xW = ■7^[Vf^{n)x{n + fi) -Vil{n- fi)x{n- fi)] 
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where is constructed using the HYP blocking trans- 
formation of Ref. . This transformation has the advan- 
tage of using only links lying within hypercubes attached 
to the original thin-link, so that is less extended than 
the fat links used in the asqtad action. HYP-blocking 
also includes SU(3) projection. We set the HYP blocking 
parameters to the values that remove the tree-level cou- 
pling of quarks to gluons having one or more components 
of momenta equal to 7r/a (with the other components 
vanishing) . In the notation of Ref. [ill, HI] , these are the 
HYP(H) parameters. These are the values we have used 
in our simulations. 

The HYP action is only partially improved — taste- 
breaking 0{a^) interactions are removed, but taste- 
conserving 0{a^) terms are not. It would thus seem to 
be a poorer choice than the asqtad action, which is fully 
O(a^) improved. It turns out, however, to be a better 
choice in practice, for two reasons. The most important 
is that it leads to substantially smaller taste-splittings 
between pions [Tsj . Since taste-splitting is the dominant 
O(a^) effect with staggered fermions, this means that 
HYP-smeared quarks have smaller 0{a^) effects than 
asqtad quarks. The second reason is that the HYP action 
is more continuum-like, in the sense that loop contribu- 
tions to matching factors are typically smaller. This is 
known explicitly at one-loop for bilinears (as found for 
the Wilson gauge action in Ref. 11] , and for the improved 
gauge action in the present work — see Table IIIip , and is 
expected also to hold at higher order because the four- 
fermion operators induced at one-loop have greatly re- 
duced coefficients compared to asqtad quarks [g^ . These 
advantages, as well as the computational simplicity of im- 
plementing this action for valence quarks, have led us to 
pursue calculations using the HYP-smeared action. For 
completeness, we note that similar reductions in taste- 
splittings (and presumably similar reductions in one and 
higher-loop contributions to matching factors) can also 
be obtained using the more highly improved, and more 
complicated, HISQ action [26| . 

The thin links are related to the gauge fields A^^ in the 
usual way, 



U^{x) = exp[igoAf_,ix + il/2)] 



(9) 



Writing the HYP links in a similar way in terms of 
"blocked gauge fields" i?^, 



V^{x) = exp 



(10) 



the Bf^ can be expressed in terms of the gauge fields as 



(11) 



n=l 



Here i?^"^ contains all terms with n powers of A. It turns 
out that we need only B^"^ in the one-loop calculation. 

(2) 

While Bli ' enters in one-loop "tadpole" diagrams, these 



contributions vanish because of the SU(3) projection 

[ll|,[23|. Thus all we need is the relationship between b'iP 
and Ay-. 



(12) 



where we have gone over to momentum space. A con- 
venient general form for the kernel h^^ (k) is (following 
Ref. but using the notation of Ref. jHI^) 



S^,^Di_,{k) + (1- 5i_,y)s^SyGu^i_i{k) ,(13a) 
l-di^sl + d2 E slsl 



(13b) 



2 7t2 



G.,^(fc) = di - d2 t """^ + + disl (13c) 



Here = sin(/c^/2). 

The coefficients di distinguish different choices of 
smeared links. 

(i) Unimproved ("thin") hnks (C/^): 

di =0, d2 = 0, c?3 = 0, d4 = 0. (14) 

(ii) HYP-smeared fat links (VJj) whose coefficients are 
chosen to remove O(a^) taste-symmetry breaking 
at tree level: 



di = 1, 



1, 



1, di = 0. 



(15) 



(iii) Fat7 links with the Lepage term (VF^), which re- 
move all O(a^) couplings at tree level, both taste- 
violating and conserving: 



di = 0, d2 ~ I, da = 1, di = 1. 



(16) 



With the exception of some tadpole diagrams, the per- 
turbative calculation requires the propagator from one 
smeared-link to another. This takes the form 

{Bl'^^\k)Bi'^'^{~k)) 
= EV(fc) h^pi-k) {Ai{k)Af,{-k)) 



= K0{k)D':i^{k) 



(17) 



where b, c are color indices, and D^^^^ is the propaga- 
tor for the improved gluon action given in eq. ([1]). In 



^ The reversal of the indices on Gv.fj, is intended and follows 
Ref. [H. 
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the third hne we have used the fact that /i^^ is an even 
function of the momentum. We note that this smeared- 
smeared propagator includes off-diagonal {fi ^ v) terms 
even if the gluon action is unimproved, because the ker- 
nel hfj^i, has off-diagonal terms. Thus, for those diagrams 
involving the smeared-smeared propagator, the general- 
ization to an improved gluon action does not introduce 
any new types of contribution, and we can carry over the 
form of most of the results from Ref. [ll| . Appendix |B] 



where y denotes the particular 2^ hypercube, and A, B 
are "hypercube vectors" denoting the positions within 
the hypercube. The matrices l/js CE5 ^f)ab are in the 
standard notation of Refs. [l^, Q^. The spin (S) and 
taste (F) of the bilinear can each be scalar, S, vector, 
V^, tensor, T^^,, axial vector, A^, or pseudoscalar, P. 
The only new feature of these operators compared to 
those used with unimproved staggered fermions lies in 
the links used to make them gauge invariant. The factor 
V'"^(y + A,y + B) is constructed by averaging over all of 
the shortest paths between y + A and y + B, and for each 
path forming the product of smeared gauge links. We use 
the same smeared links as in the fermion action, i.e. 
for the unimproved action, for the HYP action, and 
Wf^ for the asqtad action. This ensures the conservation 
of the current [V x S] for unimproved and HYP-smeared 
fermions. 

For the asqtad action, however, the presence of the 
three-link (Naik) term in the action means that [V x S] 
is not the conserved current. As a check, we have also 
calculated the matching factor for the asqtad conserved 
vector current. This is described in Appendix [Cl 

Finally, we have also implemented mean-field improve- 
ment for the HYP action ill,]. Although the HYP- 
smeared links fluctuate much less than thin links, resid- 
ual fluctuations are present and can be partly removed 
by rescaling the links. The rescaling factor, Uq^, is cho- 
sen to be the fourth-root of the plaquette constructed 
from smeared links. The details of this procedure are 
explained in Ref. [27| and we do not repeat them here. 



III. RENORMALIZATION OF BILINEAR 
OPERATORS 

The one-loop Feynman diagrams are shown in Fig. [T] 
The X and Z diagrams are infrared divergent. We regu- 
larize this divergence, following Refs. [13: El by adding a 
gluon "mass" term, A^, to the denominator of the gluon 
propagator. This allows us to set both quark masses and 



describes how this works. 

For the asqtad action, however, the situation is less 
simple, since not all the links are smeared. Some dia- 
grams must be calculated anew when using an improved 
gluon propagator. We discuss this in Appendix [Cl 

We now turn to the bilinear operators. We construct 
them from the standard hypercube convention [28^], in 
which the spin and tastes of the two continuum fermions 
are spread over a hypercube: 



(18) 

I 




(e) Z 



FIG. 1: One-loop Feynman diagrams contributing to match- 
ing factors for bilinear operators. 



external momenta to zero. 

We have undertaken two independent calculations, one 
based on the approach of Ref. [32], the other following 
Refs. [ll|,[13l- These two methods lead to identical results 
within the accuracy of the numerical integrations. We re- 
fer to these references for discussions of the methodology. 

One-loop matrix elements of lattice operators take the 



J 



[S^F]{y) = lj2iXb{y + A) {is<S>^f)ab Xc{y + B)] V'^{y + A,y + B) . 

A.B 
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general form (with the lattice spacing restored for clarity) 

^Latt,(l) ^ 
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+0(a) 



d^.-i'P log(aA) + C, 



Latt 
ij 



M 



Latt,(0) 



(19) 



where the superscript indicates the order in perturbation 
theory, and the subscript labels the different spins and 
tastes. Thus ^A/f-^'^"'^^) ig tree-level matrix element of the 

j'th bilinear operator, are the one- loop anomalous 
dimensions (which, for bilinears, are diagonal): 



= -2CFd, 



(20) 



Here Cp = 4/3, while the di depend on the spin, but not 
on the taste, of the bilinear: 

= {3,0,-1,0,3} for {S,V,T,A,P]. (21) 

Finally, C]^''** is the finite part of the correction. 
The finite part can be broken down as 



C, 



Latt 



= CpSijdi^FQQQa - 7-E + 1) 

+Cf [X^J + [Y, +T, + ZT + Z)\ . (22) 



The first line is the finite coefficient accompanying the 
log(aA), and is thus proportional to the anomalous di- 
mension matrix. The numerical values of the constants 
are Fqooo = 4.36923(1) and = 0.577216 .... The sec- 
ond line gives the finite contributions from each of the 
diagrams, and incorporates the result that only the X- 
diagrams give rise to mixing between bilinears. 

Expressions for the finite contributions are given in 
Appendices |B] (HYP fermions) andlCl fasgtad fermions). 
We present numerical values in Tables |T] (diagonal com- 
ponents) and im (off-diagonal components). Results are 
shown with both Wilson and improved gauge actions. We 
quote only two decimal places for brevity; our numerical 
evaluations are accurate to at least 1 x 10~^. For fermion 
actions, we compare the tadpole-improved staggered ac- 
tion [i.e. the action of eq. (|8]) with links U^/uo\, the 
HYP action, with and without mean-field improvement, 
and the asqtad action.^ 

We note that improving the gauge action leads to a 
moderate decrease in the magnitude of the matching 
coefficients, except for those which were already small 

{\C^^\<l). 



^ The results with the Wilson gauge action agree with those ob- 
tained in Refs. flll[27ll. with the exception of the mixing coeffi- 
cients for asqtad action [column (g) of Table |ll], where a (numer- 
ically small) error in Ref. has been found. We note that the 
off-diagonal mixing coefficients in Table llll must be multiplied by 
—3/4 to be compared to those quoted in Refs. Ilfl. l27ll . 



IV. MATCHING WITH CONTINUUM 
OPERATORS 

The continuum operators to which we wish to match 
are 



(23) 



where Q is a four-taste quark field with exact SU(4) flavor 
symmetry. 

The general form of the one-loop matrix elements of 
these continuum operators can be expressed as 



M 



Cont,(l) _ 



(4^) 



7f log(^)-^Cf™* 



M 



Cont,(0) 



(24) 



Here C'^^ont ^j^g finite part of the renormalization fac- 
tor, which, in the MS scheme using naive dimensional 
regularization for the gamma matrices, is 

Cf°"' = |y,0,|o,y| for{5,l/,r,A,P}. (25) 

Now we are ready to match the lattice and continuum 
operators. The tree-level matrix elements _A/(^^"'''°' and 
^Cont,(o) matched by construction.^ Equating the 
one-loop matrix elements in eqs. (fT9|) and (|24p leads to: 



Cont,(l) 



Latt,(l) 



(26) 



2 

+ J^[-S,,Ji"hog{^^a) + c,,] (27) 



where Zij is the matching factor at the one loop level 
and Cij is 



^Latt \ 



(28) 



A partial check of our results with the improved gauge 
action can be made by comparing to the one-loop result 
for Z„i given in Ref. [8] (as part of a two-loop calculation): 



^2 r 



1 



An 



0.1188(1) - -log(Afa) 



(29) 



For staggered fermions, one has an exact relation 

Zs = Zitgii — 1/Zm ■ (30) 



This requires that one does the unitary change of basis to convert 



from (7g igi ^f) to (75 i55 ^f) matrices, as explained in Ref. [27 
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TABLE I: Diagonal part of finite coefficients, C^i''". Indices fi, u, p and a are all different. Results are given for the following 
choices of fermion and gauge actions: (a) Tadpole-improved staggered fermions with Wilson gluon action; (6) Tadpole-improved 
staggered fermions with improved gluon action; (c) HYP fermions with Wilson gluon action; (d) HYP fermions with improved 
gluon action; (e) Mean-field improved HYP fermions with Wilson gluon action; (/) Mean-field improved HYP fermions with 
improved gluon action; {g) Asqtad fermions with Wilson gluon action; {h) Asqtad fermions with improved gluon action. For 
brevity, we quote only two decimal places of the numerical results. 



Operator 


(a) 


(b) 


(c) 


id) 


(e) 


(/) 


(ff) 


{h) 


(1® 1) 


42.47 


34.12 


3.46 


2.54 


2.06 


1.58 


6.23 


4.83 




14.86 


12.28 


0.05 


-0.24 


0.05 


-0.24 


3.77 


2.84 




2.58 


2.10 


-3.15 


-2.85 


-1.75 


-1.89 


4.40 


3.25 




-3.65 


-3.29 


-6.36 


-5.44 


-3.55 


-3.51 


6.19 


4.62 




-8.33 


-7.40 


-9.56 


-8.01 


-5.34 


-5.12 


8.37 


6.34 


ili^ ® 1) 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


-1.89 


-1.91 


(7^ ® Cm) 


6.55 


5.32 


1.57 


1.21 


0.17 


0.24 


-5.70 


-5.17 


(7m ® 


-0.23 


-0.40 


-2.46 


-1.90 


-1.06 


-0.93 


2.01 


1.32 


(7m ® 6m^) 


4.53 


3.46 


-0.49 


-0.42 


-0.48 


-0.42 


-1.91 


-2.00 


(7m ® i^p) 


-3.34 


-3.06 


-4.99 


-3.84 


-2.18 


-1.91 


5.37 


4.12 


(7m ® C-^s) 


-0.25 


-0.51 


-2.80 


-2.22 


-1.40 


-1.25 


1.37 


0.74 


(7m ® 'Jms) 


-6.52 


-5.80 


-7.58 


-5.82 


-3.36 


-2.93 


8.66 


6.87 


(7m ® 6) 


-3.69 


-3.44 


-5.29 


-4.13 


-2.48 


-2.20 


4.77 


3.58 


(7m- ® 1) 


-1.46 


-1.54 


-2.45 


-1.79 


-1.05 


-0.83 


0.72 


0.23 


(7^1- ® Cm) 


-0.42 


-0.63 


-0.50 


-0.34 


-0.50 


-0.34 


-3.79 


-3.59 


(7mi' ® Cp) 


-3.35 


-3.11 


-4.63 


-3.40 


-1.82 


-1.47 


4.90 


3.77 


(7mi' ® Cm-') 


-5.43 


-4.28 


0.94 


0.76 


-0.46 


-0.20 


-9.68 


-8.50 


(7mi' ® Cmp) 


-1.04 


-1.18 


-2.46 


-1.79 


-1.06 


-0.83 


0.82 


0.32 


(7mi' ® Cp<t) 


-5.92 


-5.26 


-6.92 


-5.10 


-2.70 


-2.21 


8.76 


7.05 



TABLE 11: 


Non-vanishing mixing 


coefficients, C,-;""". 


The notation is as in 


Table H] Note that 


mean-field 


improvement does 


not change 


off-diagonal coefficients. 












Operator-i 


Operator-j 


(a) 


(fo) 


(c)/(e) 


{d)/{f) 


{g) 


(h) 


(7m ^C--) 


(7m ® Cm) 


-4.05 


-3.33 


-0.47 


-0.43 


-1.73 


-1.49 


(7m ® Cms) 


(7m ® Ci-s) 


0.86 


0.81 


0.34 


0.33 


0.73 


0.68 


(7m ® C^fs) 


(7m 'X'Cp-'s) 


1.98 


1.72 


0.37 


0.36 


1.09 


0.98 


(7mi' ® Cms) 


(7m-' ® Cps) 


0.90 


0.73 


-0.01 


-0.004 


0.23 


0.19 



Thus Ref. (S] would predict 

Zs = l + [8 log(Ma) + C^°"' - 4.8262] , (31) 

and thus that Cl^\^ = 4.8262. This agrees with our re- 
sult, which is given in the first row of column (h) in Ta- 
ble ID 

V. DISCUSSION 

We can use our results to compare the reduction in 
the size of one-loop matching factors achieved by differ- 
ent improvement schemes. The starting point is tadpole- 
improved staggered fermions [column (a) in the Tables]. 
This comparison is most straightforward for vector (and 
axial) currents, since for these the anomalous dimensions 
vanish, so that there is no dependence on renormaliza- 
tion scale or scheme, and (7^^°"* vanishes. Thus, for these 
operators, C]^-^" gives a direct measure of the size of the 
corrections. We note that, for the lattices on which we 



are presently simulating (with a « 0.045 — 0.12 fm), the 
range of values of g^/(167r^) is 0.017 — 0.026 (evaluating 
the coupling at scale 1/a in the MS scheme).^ Thus a 
finite coefficient of size |Cj'^'^"| ~ 5 corresponds to a 10% 
one-loop correction. 

There are eight different tastes of vector currents in 
Table [H and three in Table |lll We see that corrections 
for all actions are of moderate size, with the largest mag- 
nitude being « 9, so there is not much room for improve- 
ment over the simple tadpole-improved action. We do 
note, however, that, for all fermion types, improving the 
gauge action does lead to a moderate reduction in the 
size of the correction, except when the magnitude of the 
coefficient is already of order unity. For most coefficients 
the reduction is in the range 10-25%. This reduction is, 



A similar range is obtained using a in the V-scheme evaluated 
at the scale 2/ a, which is the close to the calculated q* values 
for Zm i. 
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TABLE III: Spread of values for the diagonal finite corrections, cu, both for a given spin (leading to a scale and scheme 
dependent result), and between all operators (setting fia — 1). Notation for columns is as in Table [D 



Spin 


(a) 


(6) 




(d) 




(/) 


(3) 


(h) 


S/P 


50.8 


41.6 


13.0 


10.6 


7.4 


6.7 


4.6 


3.5 


V/A 


13.1 


11.1 


9.2 


7.0 


3.5 


3.2 


14.4 


12.0 


T 


5.5 


4.6 


7.9 


5.9 


2.2 


2.0 


18.4 


15.6 


All 


50.8 


41.6 


14.5 


12.6 


8.9 


8.7 


19.0 


16.0 



however, smaller than that achievable using mean-field 
improved HYP fermions, where the reduction is close to 
50%. Without mean- field improvement, HYP and asq- 
tad vector currents turn out to have slightly larger one- 
loop corrections than those for tadpole-improved stag- 
gered fermions. 

Turning now to the operators with anomalous dimen- 
sions, we can remove scale and scheme dependence by 
considering the differences between the cu for fixed spin 
and varying taste. We list in Table Hill the spread of val- 
ues for each of the three choices of spins. We also include 
the spread of values across all spins and tastes, choos- 
ing fi — \/a. This is a useful measure of the range of 
corrections, since the variation with ^ is relatively weak. 

We see from the table that, for HYP and asqtad 
fermions, improving the gluon action reduces all the 
spreads, although by a small amount. Once again, the 
greatest reduction is achieved by the mean-field improved 
HYP operators. We also see that for HYP operators 
without mean field improvement, the largest spread is 
somewhat smaller than that for asqtad fermions. This 
is the most important indicator when considering four- 
fermion operators, since, after Fierz transformation, bi- 
linears of all spins appear. 

Finally, it is of interest to see how other choices of 
improved gauge action compare to the Symanzik action. 
As representative examples we show, in Table ITVl results 
for diagonal coefficients for HYP-fermions with both the 
Iwasaki (c = -0.331, c' = 0) and DBW2 (c = -1.4067, 
c' = 0) [335 actions. We also repeat the results with 
Wilson and Symanzik gauge actions for comparison. We 
see that, although the changes are relatively small, both 
Iwasaki and DBW2 actions lead to smaller coefficients. 
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Appendix A: Improved gluon propagator 

The irnproved gluon propagator was determined 

in Ref. [l7[, and presented in a useful general form in 



TABLE IV: Diagonal part of finite coeflicients, C\f^^, for 
HYP fermions with (c) Wilson, (d) Symanzik, (i) DBW2 and 
(j) Iwasaki gluon actions. Other notation as in Table H] 



Operator 


(c) 


(d) 


iv 


w 


(1® 1) 


3.46 


2.54 


-1.89 


1.03 




0.05 


-0.24 


-3.07 


-0.99 




-3.15 


-2.85 


-4.18 


-2.90 


(1®Cm5) 


-6.36 


-5.44 


-5.27 


-4.78 


(1®C5) 


-9.56 


-8.01 


-6.34 


-6.62 


(7m ® 1) 


0.00 


0.00 


0.00 


0.00 


(7m ® Cm) 


1.57 


1.21 


0.44 


0.81 


(7m ® C-') 


-2.46 


-1.90 


-0.69 


-1.28 


(7m ® Cmi-) 


-0.49 


-0.42 


-0.20 


-0.32 


(7m ® i^p) 


-4.99 


-3.84 


-1.39 


-2.59 


(7m ® C-^s) 


-2.80 


-2.22 


-0.87 


-1.55 


(7m ® Cms) 


-7.58 


-5.82 


-2.10 


-3.92 


(7m®C5) 


-5.29 


-4.13 


-1.56 


-2.84 


(7m-^ ® 1) 


-2.45 


-1.79 


0.36 


-0.92 


(7m-' ® Cm) 


-0.50 


-0.34 


0.83 


0.01 


(7m-' ® Cm) 


-4.63 


-3.40 


-0.15 


-1.94 


(7m-' ® Cm-') 


0.94 


0.76 


1.22 


0.74 


(7m-' ® Cmp) 


-2.46 


-1.79 


0.36 


-0.92 


(7m-' ® Cp<t) 


-6.92 


-5.10 


-0.68 


-3.02 



Ref. [3J|. Here we present a simpler form. 

The improved gluon action takes the following 
quadratic form in the gluon fields after covariant gauge 
fixing: 

Q,^ ^ (Al) 
= -k^V^, + fk^l^ - cM^. , (A2) 
where a is the gauge-fixing parameter, 

fc"=^(fc^)" with A:^ = 2s^ = 2sin(fc^/2), (A3) 

■p is the longitudinal projector 

^M.-^ [P' = r], (A4) 

S'^ is the transverse delta-function 

C = -^M- - ['PS^ - , {S'^f = S'^] , (A5) 

M is an auxiliary transverse matrix 

A^/ii/ — S^yk^jk — k^jky — k^kj^ H r- , (A6) 
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satisfying VI^ = and 5'^ M = 0, and the function / is 

/ = (tj - c'P - ck^/k^) . (A7) 

The coefficients w, c and c' are determined from the pa- 
rameters of the improved action [see eq. ([T])]: 

w = Co + 8ci + 8c2 + 16c3 (A8a) 
c = ci - C2 - C3 (A8b) 

C' = C2 + C3 (A8c) 

I 



In the standard normahzation convention = 1, and 
this is the value we use in our perturbative calculation. 
We keep as a free parameter, however, since this al- 
lows phenomenological estimates of the impact of using 
different variants of the improved action. 

The inversion of Q is facilitated by observing that 
M.^ \s dependent on 5^ , Ai and A^^, as follows from 
the Cayley-Hamilton theorem applied to the three- 
dimensional transverse space. The result for the im- 
proved propagator is: 



P(P - CXi) +C^X2 


-t-c(fc2 - cxi)M^,y +c^{M'^)^,v 




k'^{k'^ — cxi) + '(?X2 





(A9) 



where 



tor" (|17l) . which are defined through 



xi = Tr(X) = 



(AlOa) 
(AlOb) 



X2 = [Tr2(X) -Tr(7W2)j /2 

= P [fc'^-(3/2)fc2p + (l/2)(P)3 
— 3P ^ ^ft^^^p I 



(AlOc) 



CCS = [Tr3(X) - 3Tr(>l)Tr(X2) ^ 2Tr(7W3)j /g 
{k^f 
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(fc^f -f 3(fc^)^-6fc^(fc^)^ + 8fc6fc2 



6fc* 



— 4 ^ ^ ^3 



(AlOd) 



We note that 2?^™? is symmetric, and that its off-diagonal 
elements are proportional to fcpfc,y multiplied by a func- 
tion that is even in each of the components of k. Thus it 
is convenient to write the propagator as 

= 5,.V'^^ + (1 - 5,,)k,k&;^^ . (All) 



Appendix B: One-loop results for HYP-smeared 
fermions 



In this appendix we present the one-loop expressions 
for matching factors for HYP-smeared fermions. The re- 
sults are presented in a general way such that they in- 
clude also unimproved staggered fermions, as well as the 
impact of mean- field improvement. 

Key building blocks for these results are the diagonal 
and off-diagonal parts of the "smeared-smeared propaga- 



a, 13 



= B [(5^^P,, + (1 - 5^,^)As^s^O^,y] 
Here B is the boson propagator 



B 



(Bl) 
(B2) 

(B3) 



Explicit expressions for and O^^ can be obtained 
using the decompositions (|Alip and (|13ap of the im- 
proved gluon propagator and smearing kernel, respec- 
tively. They are, however, uninformative and we do not 
reproduce them here. 

As noted in the main text, the expressions for one- 
loop matching factors for HYP-smeared fermions that 
are given in Ref. still hold as long as the and 
O^i/ defined above are used. This follows because the 
generalized and O^i, still satisfy the property of be- 
ing symmetric separately in each component of s^. This 
property is used to simplify the expressions. 

We think it useful to repeat the one-loop expressions 
here, both for the sake of clarity (since Ref. con- 
sidered other cases not relevant here), and in order to 
facilitate the subsequent discussion of the results with 
asqtad fermions. 

The diagonal part of the AT-diagram contribution is 



Xii — 



E 



clP,{s^ fBF^Vm - ^ 



(-) 



+2^ / v<s,s^O^,PFV,(fc) 
X [l - 



(B4) 
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Here the integral is 



16^2 n "^^^ 



27r 



(B5) 



The new abbreviations are Cu — cos(fcp/2) and = 
sin(fc^). For HYP fermions = s^, although this will 
not hold for asqtad fermions. The denominator of the 
fermion propagator is 



F = 



while the vertex factors are 



(B6) 



(B7) 



with and -F^ being hypercube four-vectors describing, 
respectively, the spin and taste of the bilinear (which are 
collectively labeled "«" ) . 

The non-zero mixing coefficients are (j ^ j)^ 



V / O D IT'S / N N r) -l-irm 

Xij = - J 2BF SlS2(Si S2 -P3C3VP; 



2 1 /mix 
^3 



(B8) 



S1S3 Ol2Vo.^ 



The vertex factors for the cases of non- vanishing mixing 
are collected in Table fVl 

The contribution of K-diagrams depends only on the 
distance A = '}2,^{S — i^)^. It vanishes for A = 0, and is 
otherwise 



YA = J2^k^ (A>1), 



(B9) 



fc=i 



where 



Ia= [ BF {sis^Pi + 12sls2S^02i) Vy{A) (BIO) 

Jk 



with vertex factors 



Vy{1) = 1 



Vy{2) 
Vy{3) 



C2+ C3+ C4 



C2C3 + C2C4 + C3C4 



(Bll) 



Vy(4) = C2C3C4 . 



L2 pa: 

an error in eqs. (20-23) of Ref. The numerical impact of this 
error is, however, minor: results for the mixing coefficients given 
in Table If of Ref. are changed by less than 10""^. We stress 
that for the HYP action, the expressions and numerical values 
given in Ref. Illl are correct. 



The tadpole contribution also depends only on the dis- 
tance A. It is conveniently divided into the contribu- 
tion from gluon propagators beginning and ending on the 
same smeared link, T^, and the remainder, T^, which re- 
quires A > 2. The former is naturally combined with the 
self-energy tadpole to yield 



T^ + ZT={A- 1) 



Imp - / (S/2)Pi 



(B12) 



Here Imf is present if mean-field improvement is imple- 
mented, and is given by 



^MF — 



/ Bsj [Pi 

Jk 



(B13) 



Note that for unimproved staggered fermions, mean- 
field improvement is commonly called tadpole improve- 
ment. The numerical values of Imf are tt^ = 9.869605, 
7.229736 for unimproved staggered fermions with Wilson 
and improved gauge actions, respectively, and 1.053786, 
0.722795 for HYP fermions with the same two gauge ac- 
tions. 

The second part of the tadpole contribution is unaf- 
fected by mean-field improvement, and is 

J^A = ^ 4.Bsl40i2VTiA) , (B14) 

with vertex factors Vt(0) = Vt(1) — and 
VTi2) = 1 , Vt{3) = 2 + C3 , Ft(4) = 3 + 2c3 + C3C4 . 

. 

Finally, the non-tadpole self-energy contribution can 
be obtained from the conservation of the taste-singlet 
vector current: 



Z — —Xii — Yi 



* = (7m ® 1) ■ 



(B16) 



These results hold for the HYP action with different 
choices for smearing kernel (entering through the coeffi- 
cients di-i) and different choices of gauge action (enter- 
ing through the coefficients co-3 in the gluon propaga- 
tor). 



Appendix C: One-loop results for asqtad fermions 

Using asqtad rather than HYP-smeared fermions leads 
to three changes: (i) the links are now 0{o?) improved, 
rather than HYP-smeared; (ii) the Naik term is present; 
and (iii) the hypercube vector current is no longer con- 
served. The impact of these changes is that, while the 
X- and Y-diagrams can be obtained by simple substitu- 
tions from those for HYP-smeared fermions, the tadpole 
and self-energy contributions must be calculated anew. 
In detail, the changes from the previous section are as 
follows: 

• The coefficients in the smearing kernel are now di = 
0, d2 — = = 1. These enter through and 
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TABLE V: Vertex factors for non- vanishing bilinear mixing coefficients, needed in eq. (|B8|) . The components /i, v and p are all 
different, but otherwise arbitrary. We use the shorthand = cos(fcp). Mixing coefficients for which both i and j are multiplied 
by (75 ® Cs) 8.re the same, and are not shown separately. 



Operator-i 


Operator-j 


T /mix 


T/mix 


(7m 


(7m® Cm) 


2 


2[S24^ -24'S3] 


(7m ® Cms) 


(7m ® C-s) 


-2C3C4 


-2C4[S2S^ C3 - 2s^ S3C2] 


(7m ® Cm-^s) 


(7m ® Cm-^s) 


- [C3 + C4] 


-S2S3'[C3 + C4] + 2s^ S3[C2 + C4] 


(7mi' ® Cms) 


(7m-' ® Cms) 


C4 — C3 


-S2S3' [C3 - C4] + 2S2 S-i{ci ~ C4] 



• In all expressions, now differs from due to 
the effect of the Naik term on the propagator: 



s4l + ,s2/6). 



(CI) 



The form of the result for X-diagrams remains un- 
changed, but and O^^i, are changed because of 
the impact of the Naik term on the quark-gluon 
vertex. They are replaced by P^^ and O^^, ob- 
tained from 



where 



a, 13 



(C2) 



B ['Jm.P^^^ + (1 - 5M.)45^s-,0;r/] , (C3) 



(C4) 



• The form of the result from Y-diagrams is also un- 
changed, but now and O^i, must be replaced by 



and O^j^, which are obtained from 



mp 





(C5) 



a, 13 



= B [5,,,P;^ + (1 - 5^,)A-s^-s,0^,] . (C6) 

These asymmetrical changes reflect the fact that 
the Naik term enters when the gluon attaches to 
the external fermion leg but not when it attaches 
to the operator. 

• Type-(b) tadpole diagrams are unchanged in form 
and involve O^i, without any change from the Naik 
term. In other words, one uses eq. (|B14p with O^i, 
from eq. (jB2[) . 



Type- (a) tadpole diagrams are replaced by 



T^ + Zt = (A-1) 



1 f pi-P(_5 + 3ci + 3c2 - 3ciC2) 
4 Jk 

+12 / pi-P(s-i)^(s-2)' + ^T^^ 

T'^y"'- J pi'J^%(i + ci) 



(C7) 



where the p-^i"" terms arise from the tadpole- 
improvement of the links, with 



j.Sym ^ I (2?I-P(52)2 - 4P17(S1)2(S2)' 



(C8) 



For the Wilson gauge action p-^y™ — 7^2^ ^^j^jg 
factor is reduced for the improved gauge action to 
7.229736. 



• The non-tadpole self-energy is given by 
Z = j^{B^ + BFh), 

where 

Iz = cf [1 - 2{s^ fF] [c^Pi^^ - 3c2pf ^] 



(C9) 



N uNN 



1252525^0^2"^ 



-s,s^{sl/i-cl)D[ 

-Hs?S2S^(3/4-S?)G'i,2 



(CIO) 



Here = Cp(l + s^/2), while the new quanti- 
ties P'^j and ^ arise from the propagator from a 
smeared link to a thin (Naik) link. They are defined 
by 



mp 



B 



Sf^uD' + (1 - (5p^)s^s^G' 



(Cll) 



As noted in the main text, we have also calculated the 
matching factor for the asqtad conserved vector current. 
This is given by adding 1- and 3-link terms to the hy- 
percube current, both containing tadpole-improved thin 
links: 



^xiiE''M(y) 



[V^ X S]iy;W) + -^^[V^ X S]{y;U) 



X^A^2f. Upiy + A- 2/i, y + ^+ A)Xj^A+M 



(C12) 
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where the notation for bihnears is as in eq. (I18|) . except 
that the second argument of [VJ^ x S] indicates the type of 
hnks used to create the parallel transporter. In addition, 
A is a vector running over the 8 positions of the cube 
perpendicular to fi, while "h.c." implies interchange of 
the positions of x a-nd x fields and hermitian conjugation 
of the gauge fields. 

At tree-level, and for physical external momenta, the 
extra 1- and 3-link terms in the current cancel. At one- 
loop, however, these terms lead to additional contribu- 
tions. For X-diagrams, the effect is to change the vertex 
functions as follows: 

V,{k) V-(fc) + I [V,{k) - F,(3fc)] . (C13) 



For Y-diagrams, the expression for in (IC5[) is changed 
by the substitution 

Kfi — >h^f!-S^f)C^il + c^)/2, (C14) 

(with h^^ unchanged). The T diagrams must be calcu- 
lated anew, but turn out to exactly cancel the ZT con- 
tribution (for any choice of gluon propagator). 

The net result, as we have checked analytically, is that 
the contributions of the X, Y and Z diagrams cancel ex- 
actly. This provides an important check on our result for 
the Z diagram with asqtad fermions. 
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